Properties of an algebraic spin liquid on the kagome lattice 
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In recent work, we argued that a particular algebraic spin liquid (ASL) may be the ground state 
of the S — 1/2 kagome lattice Heisenberg antiferromagnet. Furthermore, this state, which lacks a 
spin gap, is appealing in light of recent experiments on herbertsmithite (ZnCu3(OH)6Cl2). Here, 
we study the properties of this ASL in more detail, using both the low-energy effective field theory 
and Gutzwiller-projected wavefunctions of fermionic spinons. We identify the competing orders 
of the ASL, which are observables having slowly-decaying power law correlations - among them 
we find a set of magnetic orders lying at the M-points of the Brillouin zone, the familiar q — 
magnetic ordered state, the "Hastings" valence-bond solid (VBS) state, and a pattern of vector spin 
chirality ordering corresponding to one of the Dzyaloshinskii-Moriya (DM) interaction terms present 
in herbertsmithite. Identification of some of these orders requires an understanding of the quantum 
numbers of magnetic monopole operators in the ASL. We discuss the detection of the magnetic and 
VBS competing orders in experiments. While we primarily focus on a clean system without DM 
interaction, we consider the effects of small DM interaction and argue that, surprisingly, it leads to 
spontaneously broken time reversal symmetry (for DM interaction that preserves XY spin rotation 
symmetry, there is also XY magnetic order). Our analysis of the projected wavefunction provides an 
estimate of the "Fermi velocity" vf that characterizes all low-energy excitations of the ASL - this 
allows us to estimate the specific heat, which compares favorably with experiments. We also study 
the spin and bond correlations of the projected wavefunction and compare these results with those 
of the effective field theory. While the spin correlations in the effective field theory and wavefunction 
seem to match rather well (although not completely), the bond correlations are more puzzling. We 
conclude with a discussion of experiments in herbertsmithite and make several predictions. 



I. INTRODUCTION 

Recent experiments on herbertsmithite 
(ZnCu3(OH)6Cl2) have reinvigorated the longstanding 
interest in the ground state of the S = 1/2 Heisenberg 
antiferromagnet on the kagome lattice. 1-8 This material 
contains kagome layers of antifcrromagnetically coupled 
Cu 2+ S = 1/2 moments, with an exchange energy 
J as 200K. The coupling between adjacent kagome 
layers is expected to be very small. Remarkably, no sign 
of ordering - magnetic or otherwise - has been observed 
down to the lowest temperatures measured (50mK for 
some probes). Frozen magnetic moments and spin glass 
behavior are also not observed. Herbertsmithite is thus a 
candidate for the experimental realization of a quantum 
spin liquid in two dimensions. 9,10 

Various experimental probes point to a vanishing spin 
gap. Spin liquids with this property are variously (and 
equivalently) referred to as gapless spin liquids, criti- 
cal spin liquids, or long-range resonating valence bond 
(RVB) states. So far, most works studying nonmagnetic 
ground states of the kagome antiferromagnet have fo- 
cused on gapped spin liquid states, 11,12 or valence-bond 
solid (VBS) states that break lattice symmetries. 11,13,14 
A variety of microscopic calculations have provided inter- 
esting information, but are unable to establish the nature 
of the ground state. 15-24 (See also Ref. 25 for a review.) 

The experimental work on herbertsmithite led us, in 
recent work, to investigate the possibility of a gapless 



spin liquid ground state in the kagome lattice Heisenberg 
model. 26 By considering a class of Gutzwiller-projected 
fermion wavefunctions, we concluded that the variational 
ground state of the kagome model is a particular kind 
of gapless spin liquid known as an algebraic spin liquid 
(ASL). 27,28 Some of us studied the effect of a Zeeman 
magnetic field, and argued that it leads to spontaneous 
breaking of parity and XY spin rotation. 29 Gregor and 
Motrunich considered the effect of non-magnetic impuri- 
ties in the ASL, finding results consistent with the NMR 
experiments on herbertsmithite. 30 In this paper, we shall 
work out the properties of this realization of the ASL in 
more detail; this leads to a number of predictions that 
may be relevant for herbertsmithite. 

We note that Ma and Marston have considered a dif- 
ferent gapless spin liquid (Fermi surface state) using 
Gutzwiller projected wavefunctions, and have argued it 
can be stabilized by addition of further-neighbor ferro- 
magnetic exchange. 31 Also, a different route to a gapless 
spin liquid on the kagome lattice has recently been dis- 
cussed by Ryu et. al. 32 

The effective field theory describing the ASL consists 
of Nf = 4 flavors of massless, two-component Dirac 
fermions coupled to a U(l) gauge field. This descrip- 
tion is complementary to the projected wavefunction ap- 
proach - the former correctly captures universal long- 
distance properties, while the latter can provide informa- 
tion about energetics and other more microscopic prop- 
erties. Significant progress has been made recently in 
understanding algebraic spin liquids using effective field 
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theory, 27 ' 28,32-39 and it has been found that such states 
exhibit striking observable properties. The ASL is a 
quantum critical phase - like a quantum critical point, it 
supports gapless excitations and nontrivial critical expo- 
nents, but can exist as a stable zero-temperature phase 
and can be accessed with no fine tuning of parameters. 
The ASL is Lorentz invariant, and all excitations propa- 
gate with the same "Fermi velocity" vp. The symmetry 
of the ASL is much larger than that of the microscopic 
spin model, and this emergent symmetry unifies together 
a variety of superficially unrelated observables. Some of 
these observables have slowly decaying power law corre- 
lations in space and time; these are referred to as "com- 
peting orders." 

In this paper, we show how the observable properties of 
the ASL are manifested in the kagome lattice model. In 
particular, we identify the competing orders most likely 
to have slowly decaying correlations - these include mag- 
netic orders and valence-bond solid (VBS) states, as well 
as patterns of order involving vector and scalar spin chi- 
ralitics. We discuss the detection of the magnetic and 
VBS orders in experiments. We also present further re- 
sults in the projected wavefunction approach - in partic- 
ular, we give an estimate of the velocity vp, and study 
the spin-spin and bond-bond correlations of the wave- 
function. Taken together, these results inform a variety 
of predictions that may be relevant for herbertsmithitc. 
Furthermore, the estimate of vp allows us to estimate the 
low-temperature specific heat and magnetic susceptibil- 
ity, and we find reasonable agreement with experiments. 

It bears mentioning that herbertsmithitc is almost cer- 
tainly not described by a Heisenberg model alone. There 
is now significant evidence that impurities play an impor- 
tant role, especially at low temperature. 4,6 ~ 8 It has been 
suggested that antisite defects, where Zn and Cu trade 
places, constitute the dominant disorder - this leads to 
both magnetic impurities (Cu occupying Zn sites) and 
site dilution in the kagome layers (Zn occupying Cu 
sites). 6 Estimates of the concentration of magnetic im- 
purities per kagome lattice site are in the range of 5% 
to 10%. It has also been suggested that Dzyaloshinskii- 
Moriya (DM) interaction is an important perturbation 
to the Heisenberg model. 40,41 While DM interaction is 
certainly present, its magnitude is uncertain. 

Here, we have not attempted to address the effects of 
impurities and DM interaction in detail. Instead, our 
approach is to first understand the spin liquid physics in 
a clean model with only Heisenberg exchange. One can 
then include impurities and DM as perturbations to this 
idealized model. While we do discuss the effects of small 
DM interaction, the effects of impurities, and of strong 
DM interaction, are left for future work. 

We now give an outline of the paper, which also serves 
as a more detailed overview of our main results. Read- 
ers primarily interested in predictions for experiment can 
skip Sees. II and III. We begin in Sec. II A by discussing 
how the effective field theory is obtained starting from 
the Heisenberg model. The field theory can be analyzed 



in a large-A^ expansion, and our understanding of the 
physical case (Nf = 4) derives from the analysis in this 
limit. In Sec. II B, we show how the microscopic symme- 
tries are realized as transformations on the fields of the 
continuum theory. 

Section III describes our main results from the effective 
field theory approach. Section III A reviews the physics 
of competing orders in the ASL. There are two sets of 
field theory operators whose correlations are likely to 
decay slowly for Nf = 4. These arc 15 fermion bilin- 
ears, and a set of 12 magnetic monopole operators. In 
Sec. Ill B we relate each of the fermion bilinears to an or- 
dering pattern in the kagome spin model. The ordering 
patterns arising are the valence-bond solid (VBS) state 
considered before by Hastings, 13 a set of magnetic or- 
ders with crystal momenta lying at the M-points in the 
Brillouin zone, and a pattern of vector spin chirality or- 
dering, which also corresponds to one of the DM terms 
allowed by symmetry in herbertsmithite. In Sec. Ill C, 
we relate the monopole operators to ordering patterns 
in the spin model. This is less straightforward than the 
analysis for the fermion bilinears, as there is an ambigu- 
ity in determining the symmetry transformations of the 
monopoles. This ambiguity is discussed, and partially re- 
solved using algebraic relations among generators of the 
space group. 32,37-39 Because the ambiguity in monopole 
transformations can only be partially resolved, we are left 
with three free parameters describing the possible trans- 
formation laws. We make a conjecture on the value of 
two of these parameters; based on this, we find among 
the monopoles a pattern of magnetic order corresponding 
to the familiar "q = 0" state of the kagome lattice. The 
conjecture is supported by the spin correlations of the 
projected wavefunction (Sec. VI B). Depending on the 
remaining parameter, there are monopoles corresponding 
to either a pattern of spin-chirality ordering that breaks 
time reversal symmetry, or a VBS ordering pattern. 

In Sec. IV, we discuss how to detect some of the com- 
peting orders in experiment. In particular, the M-point 
and q = magnetic orders can be detected via NMR 
relaxation rate, and also by inelastic neutron scatter- 
ing - both these quantities obey universal scaling forms 
determined by the critical properties of the ASL. The 
VBS competing order can be detected via its coupling 
to an appropriate optical phonon. The lineshape of this 
phonon can be related to the VBS susceptibility of the 
ASL, which again obeys a scaling form as a function of 
frequency and temperature. 

In Sec. V, we consider the effect of a small DM interac- 
tion. Considering first DM interactions that preserve XY 
spin rotation symmetry, we argue that, surprisingly, DM 
interaction induces XY magnetic ordering, which is likely 
to be in the q = pattern. More generically, DM inter- 
action completely breaks spin rotation symmetry, and in 
this case the same argument shows that it leads to spon- 
taneous breaking of time reversal. 

In Sec. VI we turn to a further analysis of the pro- 
jected wavefunction. By construction of variational ex- 
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FIG. 1: (a) Unit cell of the kagome lattice, showing the lattice 
vectors a\ and ct2, and the numbering of the three sites within 
the unit cell, (b) Brillouin zone. The M, K and V points 
are labeled. The M points will play an important role in 
our discussion, and have thus been numbered as shown to 
distinguish among them. The reduced Brillouin zone of the 
enlarged real-space unit cell (Fig. 2) is denoted by the dashed 
rectangle. The Dirac nodes lie at ±Q, where Q — (ir/v3)y. 

cited states, in Sec. VIA we estimate the velocity v F . 
Later on in Sec. VII, this allows us to calculate the 
specific heat and magnetic susceptibility in mean-field 
theory; these quantities compare favorably with the ex- 
perimental data. In Sees. VI B and VIC we investigate 
the spin-spin and bond-bond correlators of the projected 
wavefunction. The spin-spin correlator is dominated by 
correlations in the pattern of the q = state - in real 
space these correlations fall off roughly as 1 /r 4 . Although 
the correlations of the M-point order are expected to de- 
cay more slowly than this, they are not seen; this may 
indicate that these correlations have a small coefficient 
and thus only become important at very long distances. 
The bond-bond correlation function exhibits power law 
decay, but the observed correlations are apparently domi- 
nated by Fourier components near q = and thus do not 
correspond to the Hastings VBS state (or to the VBS 
state that may arise from the monopoles, which also has 
crystal momenta at the M-points). We speculate on the 
possible meaning of this in Sec. VIC. 

We conclude in Sec. VII with a discussion of exper- 
iments on herbertsmithite, and open theoretical issues. 
Various technical details are contained in the appendices. 

II. DESCRIPTION OF ALGEBRAIC SPIN 
LIQUID 

A. Effective theory 

We are interested in the S = 1/2 Heisenberg antiferro- 
magnet on the kagome lattice, with Hamiltonian 

H = J S r ■ S r', (1) 

{rr') 

where the sum is over nearest-neighbor pairs of sites. 
The kagome lattice has a three-site unit cell; we label 
the sites by pairs (R,i), where R = n\ax + 712*12 is the 



lattice vector labeling the unit cell, and i = 0,1,2 la- 
bels the three sites (Fig. la). We choose ax = x and 
a 2 = (l/2)a3 + (\/3/2)y, so the distance between near- 
est neighbor sites is 1/2. The reciprocal lattice primi- 
tive vectors can be chosen as 61 = 2tt[x — (l/\/3)y] and 
b 2 = (47r/\/3)y, and the Brillouin zone is as shown in 
Fig. lb. 

Although our main focus is on the pure Heisen- 
berg model, we shall also consider Dzyaloshinskii-Moriya 
(DM) interactions, of the type allowed by the crystal 
symmetries in ZnCu3(OH) 6 Cl2. We refer the reader to 
Ref. 41 for the definition of the two allowed DM terms. In 
one of these, the D z -term, the DM vector points along 
the z-axis. The D p -term, on the other hand, has DM 
vectors lying in the xy-plane. 

We begin by representing the single-site Hilbert space 
in terms of S = 1/2 fermionic spinons: 

l«> = fl\0), (2) 

where a =T,|- Wc choose fermions (as opposed to 
bosons) because they allow for the description of stable, 
gapless spin liquid phases. This representation involves 
an enlargement of the Hilbert space, and must be accom- 
panied by the local constraint fl a f ra — 1 to eliminate 
unphysical empty and doubly occupied sites. 

These variables allow one to construct low-energy ef- 
fective theories for spin liquid phases, as well as corre- 
sponding trial ground state wavefunctions. The starting 
point for these constructions is a decoupling of the quar- 
tic exchange interaction using an auxiliary field residing 
on the bonds of the lattice. Neglecting the fluctuations 
of the auxiliary field, one arrives at a quadratic mean- 
field spinon Hamiltonian. In order to obtain a correct 
description of any spin liquid ground state, one needs to 
go beyond the mean-field description and consider the 
fluctuations of the auxiliary field, which play an impor- 
tant role. One way to do this is to solve the mean-field 
Hamiltonian, and then perform a Gutzwiller projection 
onto the subspace with one fermion per site; this results 
in a legitimate trial wavefunction for the spin model. Al- 
ternatively, one can recognize that the fluctuations about 
the mean-field saddle point take the form of a gauge field 
coupled to the spinons. One can then write down an 
effective gauge theory Hamiltonian, which will correctly 
capture the universal features of a given spin liquid phase. 

In a recent paper, we studied the kagome antiferro- 
magnet using Gutzwiller projected wavefunctions. 26 The 
main result of this study was that a particular spin liquid 
state, first discussed by Hastings, 13 has a very low energy, 
even without any tuning of variational parameters. This 
state has the mean-field Hamiltonian 

WMFT = E (fUr'a + H - C • ( 3 ) 

{rr') 

The hopping parameter t is sometimes written t — x J, so 
that x gives the mean-field hopping in units of J. Also, 
s rr > = ±1 encodes the background magnetic flux of the 
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FIG. 2: The doubled unit cell used to diagonalize the mean- 
field spinon Hamiltonian Eq. (3). The lattice vectors are 2ai 
and a,2, as shown. It should be stressed that the unit cell 
doubling is a gauge artifact, and does not represent a breaking 
of translation symmetry. The thick bonds are those where 
s rr — — 1, and the sites are numbered as shown. 



gauge field coupled to the spinons; it is chosen so that 
7r-flux pierces the kagome hexagons, and zero flux pierces 
the triangular plaquettes. The total number of spinons 
is chosen so that {ft a f ra ) = 1. 

Instead of a Fermi surface, Kmft has gapless Dirac 
points at the Fermi energy, near which the fcrmions obey 
a massless Dirac dispersion with velocity vf- One can 
diagonalize the Hamiltonian using the 6-site unit cell of 
Fig. 2, with the signs of s rr i chosen as shown; the Dirac 
nodes lie in the reduced Brillouin zone at positions ±Q, 
where Q = iry/\/3 (Fig. 1). One can describe the low- 
energy excitations near the nodes in terms of the La- 
grange density for Dirac fcrmions in 2 + 1-dimcnsional 
spacetime: 

£MFT = i>aa[- h^df^lpaa. (4) 

Here, rp aa is a two-component fermion field, with a =f, J, 
labeling the spin, and a = +, — labeling the two nodes 
at Q and —Q, respectively. The two components of tp aa 
correspond to the two branches of the Dirac dispersion. 
Moreover, the index fi = 0,1,2, and -f^ = (r 3 , t 2 , —t 1 ), 
where the t % are 2x2 Pauli matrices. Finally, -0 Qa = 
iip^T 3 . More details on the band structure of Hmft, as 
well as its continuum limit, and the detailed relationship 
between ip aa and the lattice spinon fields, are given in 
Appendix A. 

The important fluctuations about the mean-field sad- 
dle point specified by Hmft are encapsulated in the cou- 
pling of the spinons to a compact U(l) gauge field, and 
the spin liquid is thus an algebraic spin liquid. Because 
of this structure, this state has been referred to as the 
U(l)-Dirac state; 26 ' 29 here, however, to emphasize the 
connections with earlier work we shall refer to it as an 
algebraic spin liquid. The coupling to the gauge field is 
encoded in the following lattice gauge theory Hamilto- 



nian: 

H eS = hJ2 e2 rr'-Kj2 COS (\ VXa \p) 
(rr') P 

- fE^^/^tac.). (5) 

{rr') 

Here e rr > and the lattice electric field and vector 

potential, respectively. They reside on nearest-neighbor 
bonds of the lattice and satisfy the canonical commu- 
tation relation [a, e] — i. The electric field has integer- 
valued eigenvalues, and the vector potential's eigenvalues 
are 27r-periodic and can be taken in the interval [— n, n). 
The notation ^ p denotes a sum over all triangular and 
hexagonal kagome lattice plaquettes, and [V x a] p de- 
notes the discrete (oriented) line integral of a rr > around 
the corresponding plaquette, i.e. the magnetic flux of 
the gauge field. In general we should allow the energy 
K to differ on triangular and hexagonal plaquettes, but 
this will not be important for our discussion here. This 
Hamiltonian must be supplemented by the Gauss' law 
constraint 

(V • e) r = fUra L (6) 

where (V • e) r = X)r'~r err ' * s the lattice divergence of 
the electric field (the sum is over nearest neighbors r' of 
r). 

In the limit h^> t,K, one recovers a spin model Hamil- 
tonian. However, this effective description is most useful 
in the limit of large K, where fluctuations of the magnetic 
field are suppressed and the spinons become good vari- 
ables to describe the physics (at least for short length 
scales). The short-distance physics of the Heisenberg 
model is presumably not similar to the short-distance 
physics of H c s in the large-if limit. Instead, the idea 
is that the two Hamiltonians may have the same long- 
distance physics, i.e. they are in the same phase. 

In the large- if limit one obtains the algebraic spin liq- 
uid, which is described by the Lagrangian density 

Asff = -tpaa[- Hn(d^ + ia^lpaa 

+ ^E( £ A-ACV A ) 2 + .... (7) 

P- 

This is the so-called QED3 Lagrangian, which includes 
the minimal coupling of the gauge field to the spinons, 
as well as an explicit Maxwell term for the gauge field. 
In general we may (and must) add other perturbations 
as allowed by microscopic symmetries - such terms are 
represented by the ellipsis. While QED3 is a strongly 
coupled problem, in the sense that the interaction be- 
tween spinons and gauge field is strongly relevant in the 
RG sense, it can be understood in a large Nf limit, where 
the number of two-component fermions fields is increased 
from 4 to N f (e.g. a = 1,. . .,N f /2). The theory can 
be solved for Nf — > oo, and one can calculate pertur- 
batively in powers of l/Nf. This large- Nf expansion 
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underpins the understanding of the algebraic spin liq- 
uid fixed point, and has been discussed in great detail 
elsewhere. 28 - 33 - 35 - 36 ' 42 ' 43 

It shall be convenient to organize the four two- 
component fermions into the eight-component object 




(8) 



The t % Pauli matrices act in the two-component space of 
each Dirac fermion, so that 



r l Vt - 



(9) 



We also define a 1 Pauli matrices acting in the spin space, 
and \x % Pauli matrices acting in the "nodal" space con- 
necting the two nodes at ±Q. For example, we have 



j 3 * 




(10) 



and 



^ 3 * 



(11) 



Different types of Pauli matrices commute with one an- 
other: 



[a\^\ = [^y] = [tV] =0. 



B. Symmetries 



(12) 



For our analysis later on, we need to work out how the 
symmetries of the microscopic Hamiltonian Eq. (1) arc 
realized in the continuum theory. The procedure for ob- 
taining these results is described in Appendix B; here, 
we shall simply enumerate the symmetries and quote 
the results. The space group of the kagome lattice is 
generated by translations by the Bravais lattice vector 
R = mai + n 2 a,2 (Tr), rotations by ir/3 about the cen- 
ter of a hexagonal plaquette (R v / 3 ), and a reflection 1Z V 
mapping y — > — y. These point group symmetries are de- 
picted in Fig. 3. For defmiteness in defining the rotations 
and reflections, we take the origin of coordinates at the 
center of a hexagon. The remaining symmetries are spin 
rotations and time reversal. 

For each symmetry operation, its action on the lat- 
tice spinon fields f ra is determined by two requirements. 




FIG. 3: Depiction of the generators of the point group sym- 
metry of the kagome lattice, lZ y and R n /s. 



First, it must reproduce the correct action of the symme- 
try on the spin operator S r - this ensures that the action 
of the symmetry on all physical operators is correctly re- 
alized. This requirement does not completely specify the 
action of the symmetry on f ra , because S r is invariant 
under local SU(2) gauge transformations of the spinons - 
the symmetry may therefore be supplemented by an ar- 
bitrary SU(2) gauge transformation. However, we should 
also impose a second requirement, which is that the ef- 
fective Hamiltonian Eq. (5) be invariant. In the present 
case, this requirement will fix the symmetry operation 
completely up to multiplication by a global U(l) phase 
factor. This also determines the transformation proper- 
ties of the electric field e rr ' and vector potential a rr >. 
Due to the presence of the gauge structure, symmetries 
have nontrivial action on the spinon fields f ra ; the math- 
ematical structure describing the realization of symme- 
tries in such a situation is referred to as the projective 
symmetry group (PSG). 44 

If S is a space group symmetry mapping r — ► S(r), 
then the above requirements dictate that the fermions 
transform as 



-fs(r)c 



(13) 



where 7r r is a gauge transformation depending on S; in 
the gauge chosen in Eq. (5) we can take n r = ±1. Spin 
rotations send f ra — ► U a pf r p, where U is an SU(2) ma- 
trix. Time reversal acts as follows: 



r-.fr 



(io- 2 )af3fr/3- 



(14) 



These transformation laws also imply that the electric 
field and vector potential transform as scalars under 
translations, and as vectors under rotations and reflec- 
tions. Under time reversal, the electric field is even and 
the vector potential is odd. 

Following the procedure in Appendix B, we obtain the 
following transformations for the (real space) continuum 
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field ^(r): point r. The meaning of Eq. (22) is that the quantities on 

the left and right hand sides have the same correlations, 

T ai : * - (V)* (15) 

at distances much greater than a short-distance cutoff. 
7a 2 : * — > (ip, 3 )^ (16) For any operator in the field theory, generically the co- 

_ 1 /%Tx \ efficient will be nonzero if and only if d transforms 

Tty ■ ^ — > (ir ) cxp ^— ^, ry J ^ (17) identically to m(r) under microscopic symmetries (in this 

/i7r \ / 27rz \ case, spin rotations, time reversal and space group oper- 

R-it/3 '■ * — » ex P ("g" T J CX P ("a - ^/* ations). Apart from this condition, the a are nonunivcr- 

_. _ / ?w - ?w • ?\ T / -i r\\ sal. This tells us, then, that the dominant long-distance 

v /v /v ^ ' v y correlations of m(r) are given by the d appearing m 

w here Eq. (22) with the smallest scaling dimension A; so, for 

^ example, (m(r)m(O)) cx |r|~ 2A . 

jj, — (p, 1 + fi 3 ) (20) Which operators in the ASL critical theory have small- 

est scaling dimension? In the large- Nf limit, these are 
1/12 3\ (91 \ a set of fermion bilinears transforming as the adjoint of 

^ R ~ 7^ M ~ SU(iV». In the physical case of AT/ = 4, these operators 

are 



III. MANY COMPETING ORDERS: FERMION 
BILINEARS AND MAGNETIC MONOPOLES 

A. General discussion 

The focus of this paper is on observable properties of 
the algebraic spin liquid on the kagome lattice; the most 
dramatic such property is the unification of many com- 
peting orders within the ASL. 36 This is manifested in 
the fact that correlations of a variety of superficially un- 
related observables obey power law decay in space and 
time with the same critical exponent. This unification 
is accomplished mathematically through the presence of 
an emergent SU(4) symmetry that contains the micro- 
scopic SU(2) spin rotation symmetry as a subgroup. This 
means, for example, that magnetic order parameters can 
be rotated into nonmagnetic ones. Furthermore, because 
the ASL is an interacting critical state, the critical expo- 
nents governing many physically interesting correlation 
functions are not equal to their mean field values. In 
particular, for those observables we refer to as compet- 
ing orders, the correlation functions decay more slowly 
than in mean field theory. Such slow decay of a correla- 
tion function indicates that the system is closer to being 
ordered in a particular channel, and is likely to be ob- 
servable in both experiments and numerical studies. 

In the language of critical phenomena, the dominant 
long-distance correlations within the ASL, and hence the 
dominant "competing order parameters," are given by 
finding the operators in the critical theory with smallest 
scaling dimension. Suppose we are interested in the cor- 
relations of some microscopic observable m(r), which is 
some function of spin operators S r > for r' near the lattice 
point r. Its long distance correlations can be understood 
from the following expression: 

mW-^cAW. (22) 

i 

On the right hand side of this expression are a set of op- 
erators in the field theory d, also located at the spatial 



N a = & T 3 T a y, (23) 

where a = 1, . . . , 15 and the T a are the 15 generators of 
SU(4). One can choose a basis for the generators where 
they are expressed in terms of the // and a 1 Pauli matri- 
ces: 

T° = {aV,<7V}- (24) 

In the large- Nf limit, Rantner and Wen calculated the 
scaling dimension of one of the N a in the context of the 
"staggered flux" spin liquid on the square lattice. 28 In 
Ref. 36 the SU(4) symmetry was emphasized, which im- 
plies that all the N a have the same scaling dimension. 
The result of Ref. 28 is 

32 

A N ^dimN a = 2-^ w -+0(l/N]). (25) 

Because the U(l) gauge interaction, which is stronger for 
smaller Nf, tends to bind together the \E" and ^ fermions 
in N a , and make the mode it creates propagate more 
coherently as opposed to decaying into its two constituent 
fermions, it is expected that Ajy < 2 for all values of Nf. 

There are other operators in the ASL that may have 
even smaller scaling dimension than the N a when Nf = 
4. The most likely candidates are magnetic monopolc op- 
erators - these are topological disorder operators for the 
U(l) gauge field. Because the U(l) gauge field is com- 
pact, its magnetic flux is quantized in multiples of 27r, 
and monopole operators are those which insert integer 
multiples of 27r flux. Therefore they carry nonzero inte- 
ger charge under the U(l)fl ux symmetry, which emerges 
at the ASL fixed point and corresponds to the irrele- 
vance of monopole operators there. 35 It is important to 
emphasize that, for the ASL fixed point to be stable, all 
monopoles that are allowed as perturbations to the ac- 
tion by microscopic symmetries must be irrelevant. The 
monopoles we will consider here carry nontrivial quan- 
tum numbers and are not allowed perturbations, so they 
may become relevant without destabilizing the ASL. 
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As is typical for topological disorder operators, it is 
difficult to construct monopole operators explicitly in 
terms of the fermions and gauge field. However, they 
can be constructed by exploiting the state-operator cor- 
respondence from conformal field theory, and it has been 
found 33 that monopoles with smallest scaling dimension 
at large-TV/ have unit magnetic charge and transform 
under the completely antisymmetric self-conjugate rep- 
resentation of the SU (Nf) flavor symmetry. In the case 
Nf = 4 this representation is 6 dimensional. Formally 
we can represent these six operators by m*^ = — m*^ a , 
where a, (3 — 1, . . . , 4 are SU(4) indices. There is another 
set of six monopoles with charge —I and the same scal- 
ing dimension, m a p = (m^)^. In the large- iV/ limit the 
scaling dimension of these operators was found to be 

A m = dim m aj8 = dV) + 0(1), (26) 

where c w 0.265. 33 If we naively put Nf — 4 into the 
leading order large- Nf result, we obtain A m w 1.06, com- 
pared to A N w 1.73 putting N f = 4 into Eq. (25). So it 
may well be the case that the monopoles give the dom- 
inant long distance correlations for the physical case of 
N f =4. 

Below, we shall work out the quantum numbers of N a 
and m Q( 3 under microscopic symmetries, and determine 
the spin model observables to which they correspond. 

B. Fermion bilinears 

It is convenient to break the N a into three different 
classes of operators. These are 

N l A = I'VVer* (27) 
N B = fVo-f (28) 
N l c = ^tV*. (29) 

Using the results of Sec. II B we can determine how each 
of these observables transforms under all microscopic 
symmetries. Clearly 7V^ form a set of 3 spin triplets, 
and Nb is also a spin triplet, while the N c are spin sin- 
glets. Under time reversal N\ is odd, while Nb and N c 
are even. This implies that N\ is the order parameter 
for a magnetically ordered state. 

For each of these operators, we can identify one or more 
microscopic operators that transform identically under 
all the symmetries of the kagome model. Then, by Eq. 22, 
the N a will contribute to the long-distance correlations of 
these microscopic operators. Actually, rather than look- 
ing for operators, it is easier to look for ordering patterns 
(i.e. states as opposed to operators) with the correspond- 
ing symmetry transformations. To understand how to do 
this, let us consider the transformations of N\. If S is a 
space group operation, then we have 

S:N\^U£(S)N{, (30) 

where U Fl (S) is an 0(3) matrix. These matrices form 
the Fi irreducible representation of the space group (see 



Appendix C for more details). If we can find a set of 
magnetically ordered states for which (S r ) transforms 
under the same representation of the space group, then 
N\ is an order parameter for this state. Furthermore, we 
can go on to construct microscopic operators that are also 
order parameters for this state, and identify correlation 
functions of these operators that obey power law decay 
with exponent 2 An. 

The relevant details of the group theory and repre- 
sentation theory of the kagome space group are summa- 
rized in Appendix C. As stated above, each component 
in spin space of N\ transforms in the F± irreducible 
representation. The N c also transform as F\, that is 
S : N' c — > U^ 1 (S)N C , and each component in spin space 
of Nb transforms as the Ai representation. 

We have already established that the N\ correspond to 
magnetically ordered states. We are primarily interested 
in finding a combination of spin operators corresponding 
to each N A , and for this purpose it is enough to consider 
collincar ordering patterns only. We therefore focus on 
the ^-component in spin space (N A ) Z . 

Our task is then to find three ordering patterns of 
kagome spins pointing along the z-axis in spin space. 
One approach is simply to guess the right pattern, but it 
is possible to be more systematic. We are interested in 
patterns of site ordering on the kagome lattice that are 
invariant under translations by 2a\ and 2a 2 , since the 
N a are also invariant under such translations. So we can 
consider all possible site ordering patterns on a cluster 
of 2 x 2 unit cells with 12 sites and periodic boundary 
conditions. For each site of the cluster r we can asso- 
ciate a vector |r), and the various site ordering patterns 
can be represented as real linear combinations of these 
states, where the coefficients of \r) should be associated 
with (S*). The space group acts on this vector space 
by S\r) = |5(t")), where S is a space group operation, 
so the vector space is a 12 dimensional representation 
of the space group. This can be decomposed into irre- 
ducible representations using the character table of Ap- 
pendix C. The irreducible representations we find in this 
manner are a complete catalog of kagome site ordering 
patterns invariant under T 2ai and T 2a2 . We find that the 
Fi representation of interest occurs precisely once in this 
catalog. The corresponding magnetic ordering patterns 
can be read off from the basis vectors of this representa- 
tion, and are shown in Fig. 4. We stress that this is the 
unique site-ordering pattern corresponding to (N A ) Z . 

Next, we wish to find ordering patterns corresponding 
to Nq. Because these operators are spin singlets and 
invariant under time reversal, it is natural to look for 
bond ordering patterns, with order parameters that can 
be built from — (S r ■ S r >), where r and r' are nearest 
neighbors. This observable measures the strength of sin- 
glet formation on a bond. Following the same procedure 
as for the site ordering patterns, we note that the same 
unit cell contains 24 bonds; the resulting 24-dimensional 
vector space is decomposed into irreducible representa- 
tions in Appendix C. In this case, we find the F\ repre- 



FIG. 4: Magnetic ordering patterns corresponding to (N\) z . 
These can be thought of as stripes of up- and down-spin 
hexagons, which are labeled with "+" and respectively. 
A spin between two up-spin hexagons points up, and between 
two down-spin hexagons points down. Spins between an up- 
spin and down-spin hexagon have no average moment. The 
three patterns are labeled according to their crystal momenta, 
which lies at the three M-points of the Brillouin zone Mi, 
where i = 1, 2, 3. The labeling of the M-points is as in Fig. 1. 

sentation occurs twice, and leads to the patterns shown in 
Figs. 17, 18 and 19. Taken together with the "uniform" 
state, where (S r ■ S r '} is the same on every bond, we can 
superpose the two patterns of Fig. 19 to form the Hast- 
ings VBS state, as shown in Fig. 5. It has three inequiva- 
lent bonds in its unit cell, and this is precisely because it 
is built from a superposition of three bond ordering pat- 
terns belonging to distinct irreducible representations of 
the space group. The fourfold degenerate Hastings state 
is not the only ground state that can be built from the 
bond ordering patterns obtained in Appendix C, but we 
restrict our attention to it for simplicity. Furthermore, 
the ordering patterns obtained in Appendix C contain 
enough information to work out the contributions of Nq 
to the bond-bond correlation function within the ASL. 

Finally we turn to Nb- This object is a spin triplet, 
but is even under time reversal; the simplest micro- 
scopic operator with these properties is the vector chiral- 
ity C rr ' = S r x S r > defined on nearest-neighbor bonds. 
Rather than searching for states as above, in this case 
it is simpler to find a spin operator transforming identi- 
cally to Nb- We define an object that naturally resides 
on the hexagonal plaquettes of the kagome lattice, which 
we label with the position vectors tv 

C h (r h )= ^ Crr,. (31) 

{rr')£h 

This can be interpreted as the vector spin chirality of a 
kagome hexagon. Here, the sum is over the 6 bonds con- 
tained in the perimeter of the hexagon, following the ori- 
entation convention shown in Fig. 6. More precisely, r is 



FIG. 5: Depiction of the Hastings valence bond solid state. 
There are three inequivalent sets of bonds, these are distin- 
guished by shading and vertical or horizontal hatching. In 
previous treatments in the literature, 13,26 the shaded bonds 
were taken to have strong dimerization, and the vertically and 
horizontally hatched bonds to have weak (and equal) dimer- 
ization. However, other choices are possible without changing 
the symmetry of the ground state. 




FIG. 6: Orientation of kagome lattice bonds used to define 
the vector spin chirality of a hexagon, Ch(r fe ). 



always taken to be the the "first" site on the bond accord- 
ing to Fig. 6, and r' the "second" site. It is straightfor- 
ward to check that Ch{i*h) has identical transformation 
properties as Nb (once a suitable long-wavelength aver- 
age is taken). The observable Cn is also related to the 
Dzyaloshinskii-Moriya (DM) interaction on the kagome 
lattice, and, in particular, 

HvM = D z Y,C z h {r h ) (32) 

r h 

is precisely the D z -term allowed by the crystal symme- 
tries in ZnCu 3 (OH) 6 Cl 2 41 . The effects of DM interaction 
on the ASL are discussed further in Sec. V. 
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C. Magnetic monopoles 

Here, we work out the quantum numbers of the mag- 
netic monopole ™*«, which was formally defined in 
Sec. Ill A, based on the results on Ref. 33. Our analysis 
builds on results of Refs. 32,37-39, where monopole quan- 
tum numbers were worked out in different but closely re- 
lated settings to the present one. However, we adopt a 
different perspective; we believe this clarifies some of the 
issues involved, and we comment on this at the end of 
this section. Our strategy is first to determine how each 
microscopic symmetry is embedded in the full symmetry 
group of the low-energy effective theory. We do not have 
enough information to determine this completely, so sev- 
eral free parameters have to be introduced. To proceed 
further, we argue that these free parameters must sat- 
isfy certain constraints, determined by relations among 
the generators of the space group. After taking the con- 
straints into account, we shall be left with three free pa- 
rameters, and we make a conjecture that determines two 
of them, based on a physical argument and on calcula- 
tions using the projected wavefunction. 

First, it is useful to recall how m* a transforms under 
the continuous symmetries of the low-energy ASL fixed 
point. It is a scalar under Lorentz and continuous trans- 
lation symmetry. Under SU(4) rotations we have 



m * a — > U ai Ufism* lS , 
and under a U(1)a ux rotation we have 



m af 3 



e m a(3 

e~ ie m af} . 



(33) 

(34) 
(35) 



We now decompose the SU(4) symmetry into its 
SU(2) x SU(2) subgroup. The first SU(2) is simply spin 
rotation symmetry (generated by the a 1 ), and the sec- 
ond consists of "nodal" rotations (generated by the fi 1 ) 
that mix the two Dirac nodes but commute with spin 
rotations. This decomposition is useful because it sepa- 
rates spin rotations from space group operations, which 
are realized in the SU(4) space as nodal SU(2) rotations. 
We replace the SU(4) index a by a pair of SU(2) indices: 
a — > (era). Here, a = 1,2 transforms under spin SU(2), 
and a = 1, 2 under nodal SU(2). Under an SU(2) x SU(2) 
rotation, we have 



l (aa)(r,b) 



where II s and U N are SU(2) matrices in the spin and 
nodal spaces, respectively. We can also decompose the 
monopole operators according to their transformations 
under the SU(2) x SU(2) subgroup. We have 

v* = [(i(T 2 )c7 l } an (ifi 2 ) ab m* {aa)(vb) (37) 
w* = (icr 2 )^[(iM 2 )^]a6m^ a)(?)f)) . (38) 

Here, v* is a spin triplet and a nodal singlet, and w* is a 
spin singlet and nodal triplet. 



We now wish to specify the action of microscopic sym- 
metries on v* and w* . In doing this, we encounter two 
problems. First, we determined the SU(4) rotation corre- 
sponding to each microscopic symmetry by working out 
the transformations of the fermion field 'P. This only de- 
termines the SU(4) rotation up to multiplication by the 
matrix C4 = diag(i, i, i, i), which generates the Z 4 center 
of SU(4). The reason for this ambiguity is that multi- 
plication of 'J by C4 is indistinguishable from a gauge 
transformation. However, m a is odd under C4, but is of 
course gauge invariant. Second, each microscopic sym- 
metry may come along with a rotation in the U(l)a ux 
space. Because \& carries no U(1)a ux charge, we have 
no information about this rotation. Both of these am- 
biguities can be taken into account by multiplying m a p 
by an undetermined phase factor associated with each 
symmetry operation. We note that there are no such 
ambiguities associated with continuous spin rotations - 
it can be shown that adding an additional U(l) phase fac- 
tor to a spin rotation is not consistent with spin rotation 
symmetry. 

We can now enumerate how the space group symme- 
tries act on the monopoles: 



T 


: v l - 




(39) 


T ai 


w l - 


- e^i%(T a >* 


(40) 


TZy 


: v* - 


-> e^Vi 


(41) 


TZy 


< - 


- e i ^R ij {n y )w j 


(42) 


R-k/3 




■* e^v* 


(43) 


Rtt/3 


Wi - 


- e^R l3 {R,^)w*. 


(44) 



Note that the reflection symmetry sends v* to Vi, and 
similarly for w* . This is because reflections change the 
sign of the magnetic charge. Here we have introduced 
SO(3) matrices describing the rotations in the nodal 
SU(2) space, which are obtained from the symmetry op- 
erations of Sec. II B and are 



R(T ai ) = 

R(ll y ) = 

R{R-K/i) = 



-1 











1 





















- 


-1 


i) 


1 










1 











-1 


-1 









(45) 
(46) 
(47) 



We can immediately eliminate the phase factor <fi ry by 
redefining v* — > e l< ^ ry ^ 2 v*, and similarly for w* . This 
does not affect the other symmetries, and we have 



TZy : v* 
1Z„ : w* 



R lj (TZ v )w j . 



(48) 
(49) 



We shall now determine the unknown phase factors 
<p ai and 4>r to the greatest extent possible, exploiting re- 
lations among the space group generators. For example, 
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we shall demand the relation R e ,„ = 1 holds when acting: 

TT j 6 ° 

on monopole operators. In general, such relations need 
only hold up to a phase when acting on quantum states. 
However, in the present case, these relations hold with 
no additional phase factors when acting on any physi- 
cal state in either the spin model or the effective lattice 
gauge theory Hilbert spaces. Since nothing forbids the 
insertion of a single monopole in the lattice gauge theory, 
the relations must hold for the monopole operators with 
no extra phase factors. The relation R%j 3 = 1 implies 
that (f>R — 7rriij/3, where tir is an integer. Next, the 
relation 

Rw /^T^Rn /?,T~^ R V J 3 = T ai (50) 

implies <fr ai = 0. So we are left only with the free param- 
eter n R — 0, 1, 2, 3, 4, 5. 

We also need to consider time reversal. We have 

T : v* — > t v Vi (51) 
T : w* -> t w Wi, (52) 

where t v ,t w = ±1. Note that time reversal changes the 
sign of magnetic charge, so it takes monopoles to anti- 
monopolcs. Furthermore, time reversal commutes with 
the SU(2) x SU(2) subgroup of SU(4) [see Eq. (19)], so 
it must take the spin triplet v* to the spin triplet Vi, 
and similarly for the nodal triplets w* and Wi. It is not 
consistent for t v and t w to be arbitrary phases, given 
that T 2 = 1, which holds for all physical states of the 
spin model and of the effective lattice gauge theory. The 
relations of Eq. (51) and Eq. (52) are the most general 
transformations consistent with these considerations. 

We proceed by using some physical input to conjecture 
the likely values of and t v - the resulting conjecture 
is supported by the numerical results on the projected 
wavefunction (see Sec. VI B). The source of physical in- 
put is the following puzzle about the kagome ASL: in 
the large- S limit of the kagome antiferromagnet, there 
are an infinite number of nearly degenerate magnetically 
ordered ground states. These consist of all states where 
the vector sum of classical spins on every triangle is zero 
- among these states, those with coplanar spin configu- 
rations are selected by the zero point energy of spin wave 
fluctuations about the ground state. Finally, among the 
coplanar states, the y/3 x \/3 state is selected by anhar- 
monic fluctuations. 45-47 Up to this point in our analysis, 
we have not found signatures of any of these classical 
low energy states within the ASL. (The ordered states of 
Fig. 4 have distinct symmetry from the classical ground 
states lying at the M points, which transform in the F 3 
representation of the space group.) It is possible that no 
sign of the large- S physics survives down to S = 1/2, but 
it would perhaps be more natural if at least a hint of it 
remained, especially given that the ASL has no spin gap. 

This leads us to the conjecture that Vi will be the or- 
der parameter for one of these classical ground states. We 
know from the analysis above that carries zero crys- 
tal momentum, and the unique classical low-energy state 




FIG. 7: The q = magnetically ordered state. 



with this property is the so-called q = state, shown in 
Fig. 7. It is the classical ground state for ferromagnetic 
third-neighbor exchange J 3 < 0. 45 It turns out that if we 
choose hr — 2 and t v = —1, then Vi is an order parame- 
ter for the q = state. This can be seen by following the 
construction of the q = state order parameter in Ap- 
pendix D. Based on these considerations, from now on 
we shall fix tir — 2 and t v = — 1. This choice is supported 
by the presence of substantial q = spin correlations in 
the projected wavefunction. 

The parameter t w = ±1 still remains to be deter- 
mined. The Hcrmitian operators wf — Wi + w* and 
= i(u>i — w*) are spin singlets, with crystal momenta 
lying at the M-points. Ift w = l, they are even under time 
reversal, and then likely correspond to bond ordering pat- 
terns. If t w = — 1, they are odd under time reversal, and 
correspond to ordering patterns in the scalar spin chiral- 
ity Si ■ (S2 x S3), where the spins can be taken on any 
three distinct nearby lattice sites. Since the bond-bond 
correlations in the projected wavefunction are apparently 
dominated by Fourier components near q = 0, we are led 
to speculate that t w = — 1 and the monopoles correspond 
to spin chirality order. 

We now contrast the approach taken here with that 
of Refs. 32,37-39, where symmetry transformations of 
monopole operators were obtained by construction of 
mean-field ground state wavefunctions with a back- 
ground ±27r flux - this is the quantum state one ob- 
tains after insertion of a monopole operator. The cor- 
respondence of the resulting states to scaling operators 
of the low-energy critical theory, in which one is ulti- 
mately interested, is not clear, and this issue makes it 
difficult to interpret the results of Refs. 32,37-39. How- 
ever, our approach here, which focuses directly on the 
scaling operator m* af} , is formally equivalent to the anal- 
ysis of Refs. 32,37-39, and puts it on firm ground. On the 
other hand, Refs. 32,37-39 obtained additional informa- 
tion on symmetries by carrying out an adiabatic insertion 
of 2ir flux numerically. In our opinion, because the ASL 
has gapless excitations, the correspondence between this 
procedure and insertion of the scaling operator to*^ is 
not established. 
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IV. DETECTION OF COMPETING ORDERS 

The competing orders arising in the kagome ASL can 
be detected by a variety of experimental probes. Here, 
we outline some predictions that we hope will be tested 
in future experiments on ZnCu3(OH) 6 Cl2. In this sec- 
tion we focus on the case of an ideal system free of per- 
turbations such as impurities and Dzyaloshinskii-Moriya 
interaction. We postpone discussion of these issues to 
the following section, and also Sec. VII. 

Two critical exponents enter into this discussion. The 
first is rjN = 2Ajv — 1, which characterizes the correlations 
of the N a fermion bilinears, and hence the M-point mag- 
netic order, the Hastings VBS order, and the hexagon 
vector chirality Ch{rh)- The second, r\ m — 2A m — 1, 
characterizes the m a p monopole operators, and hence the 
correlations of q = magnetic order. 

The magnetic competing orders can be detected via 
neutron scattering and NMR. In neutron scattering, the 
structure factor will exhibit critical scaling behavior near 
the T and M points in reciprocal space, with the scaling 
form 

*« (9 + g„) = ^(^,^), (53 , 

where c is a nonuniversal constant prefactor, and \q\ is 
much smaller than the Brillouin zone size. r\ = t]n if Q 
lies at one of the M points, and rj = r\ m if Q = 0. There 
are two different universal scaling functions Fq(x,j/), 
again depending on whether Q lies at the T point or 
one of the M points. 

Based on Eq. (53), the NMR relaxation rate is pre- 
dicted to be 

-LcxT". (54) 

Here, rj depends on the nuclear site considered - Cu, O 
and CI sites are all sensitive to the M-point magnetic 
order, but only Cu and O are sensitive to q = order. 
Therefore we expect 

VGu = rjo =mm(ri N ,ri m ) (55) 
Vc\ = m- (56) 

The VBS competing order can be detected via its cou- 
pling to phonons. We consider, for simplicity, a single op- 
tical phonon that couples to one of the patterns of VBS 
order shown in Figs. 17 and 18, and hence to some linear 
combination of the Nq fermion bilinears of the ASL. Us- 
ing inelastic X-ray or neutron scattering, the lineshape 
of this phonon can be measured. The lineshape is de- 
termined by D(uj), the retarded Green's function for the 
phonon mode's normal coordinate. We treat the phonon 
mode using the RPA of Ref. 48, which neglects the in- 
fluence of the gapped optical phonon on the gapless spin 
system, and has been successfully applied to quasi-one- 
dimensional magnets, for temperatures above the spin- 



Peierls transition. 49 We then have, for the phonon spec- 
tral function A(ui) = — (1/n) Imfl(w). 

= -g 2 XvBs(^^ r )/ 7r 

P - fig - <? 2 xVbs(^ T)} 2 + [ 5 2 Xvbs(^ T)} 2 ■ 

(57) 

Here, g characterizes the spin-phonon coupling, ilo is the 
bare phonon frequency, and Xvbs an d Xvbs are the- real 
and imaginary parts, respectively, of the retarded Green's 
function xvbs of the VBS order parameter. We have the 
scaling form 

Xvbs(-,T) = ^F VB s(^), (58) 

where c is a nonuniversal prefactor. Also, F V bs is re- 
lated to the scaling function appearing in Eq. (53) by 
lmFyBs( x ) — Fq( x ,0), for Q lying at one of the M 
points. It may be possible to test this scaling form by 
measuring the T-dependence of A(u>). 



V. DZYALOSHINSKII-MORIYA INTERACTION 

We now consider the effect of a small DM interaction, 
D z ,D p <C J- First, we must understand which new per- 
turbations to the ASL fixed point are allowed by the 
now reduced microscopic symmetry. While spin rotation 
symmetry is completely broken, the space group sym- 
metry remains unchanged. However, lattice reflections 
and rotations must now be accompanied by appropriate 
operations in spin space. In the three-dimensional Her- 
bertsmithite structure, the reflection lZ y is realized as a 
7r-rotation about the x-axis (passing through the center 
of a hexagon), along with a corresponding rotation in 
spin space. The rotation R v / 3 is realized by first mak- 
ing the same 7r-rotation about the x-axis, followed by 
mirror reflection through the plane with normal et2 — a\ 
(intersecting the center of the same hexagon). In the 
continuum, the resulting modified symmetry operations 
are 

n' y : * -» (ur 1 )(tT 1 )e < ' r ' , -»/ 2 tf (59) 

R l . y _^ e i™ 3 /3 e «rT 3 /6 e 27rW3^ (g ) 

Using these operations, combined with translations 
and time reversal (which arc unchanged) , it can be shown 
that the only fermion bilinear allowed by symmetry is 
Ng. Furthermore, the q = 1 monopole operators m*^ are 
still forbidden by symmetry - the spin triplet monopolcs 
Vi are odd under time reversal, and the spin singlet 
monopolcs Wi carry nonzero crystal momentum. There- 
fore we expect iVJ, to be the most relevant perturbation 
to the ASL generated by the DM terms. 

We now consider adding the term 

£ DM = mN z B (61) 
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to the ASL fixed point. We wish to understand how m 
depends on D z and D p . As noted in Sec. IIIB, the D z - 
term and Ng are symmetry-equivalent, so m will con- 
tain a term linear in D z . We now give an argument 
that m contains no term linear in D p . Consider a mi- 
croscopic Hamiltonian with given values of D p and D z . 
Upon some amount of coarse-graining, we obtain a con- 
tinuum field theory with a given value of m. Now we 
make a 7r-rotation in spin space about the z-axis. This 
sends D p — > — D p , but otherwise leaves the microscopic 
Hamiltonian unchanged. Various operators in the contin- 
uum field theory will be affected by this operation, but 
the value of m remains the same. We have shown that 
m does not depend on the sign of D p . Therefore, up to 
second order in D z and D v , 



D 2 

m = c zl D z + Cz2—j- 



D 2 

v 

J ' 



(62) 



where c z i, c z2 and c p2 are dimcnsionless coefficients. 

Ignoring coupling between the gauge field and the 
fcrmions, the effect of the perturbation Eq. (61) is to open 
a spin gap. However, coupling to the gauge field plays 
a surprising and important role. As the fermions are 
now massive, magnetic monopoles will condense, leading 
to confinement. 50 Naively this leads to a fully gapped 
spectrum, but we argue below that, in this case it leads 
instead to XY magnetic order. 

Let us first imagine D z ^ 0, but D p = 0, so that we 
have U(l) spin rotation symmetry about the z-axis in 
spin space. The mass term £dm will be induced, and we 
observe that this term can be written 



-DM 



m 



(63) 



That is, this term has a mass m for the up-spin fermions, 
and a mass of opposite sign — m for the down-spin 
fermions. Ignoring coupling to the gauge field for the 
moment, such a mass term leads to a v = 1 quantum 
Hall effect (QHE) for the up-spin fermions, and a v = — 1 
QHE for spin-down fermions. 51 

Now, imagine we adiabatically insert a localized +2n 
flux of the gauge field. The QHE implies that a sin- 
gle extra spin-up spinon will be induced along with the 
gauge flux, while one spin-down spinon will be depleted. 
This is equivalent to a spin flip operation, and we see 
that insertion of a monopole is accompanied by a spin 
flip. Furthermore, because the fermions are gapped, the 
gauge field dynamics is controlled by the Maxwell term 
^Maxwell = ^2 (^x9 v a\) 2 , and therefore insertion of 
a monopole only costs finite action in the imaginary time 
functional integral. This implies that the monopole prop- 
agator is long-ranged - that is, if m*{r) is a monopole 
creation operator, then (m*(r)m(0)) approaches a con- 
stant as \r\ — > 00. Because monopole insertion is accom- 
panied by a spin flip, in the ground state we will have 
(S + ) 7^ 0, which corresponds to XY magnetic ordering. 
Furthermore, within the resulting ordered state, the pho- 
ton of the gauge field is expected to correspond to the 



XY Goldstone mode. We note that a similar situation 
arises, in the absence of DM interaction, when a Zeeman 
magnetic field is applied. 29 

In Sec. IIIC, we argued that some of the monopole 
operators in the ASL correspond to the q = ordered 
state. We therefore expect that the XY order induced 
in the presence of the D z term is in the q = pattern, 
and therefore that the ASL is unstable toward q = 
magnetic ordering in the presence of D z . We remark 
that, once D p ^ 0, the U(l) spin rotation symmetry is 
lost and there is no longer an XY Goldstone mode. On 
the other hand, time-reversal is still a good symmetry of 
the Hamiltonian, but is broken in the q = state. So, 
for more general DM interaction (D p ^ 0), we expect the 
ASL is unstable to a time-reversal broken ground state. 



VI. ANALYSIS OF PROJECTED 
WAVEFUNCTION 



We now turn to an analysis of the Gutzwillcr pro- 
jected variational wavefunction for the ASL. To obtain 
this wavefunction, one begins with IV'mft)' the mean ~ 
field ground state of the Hamiltonian Eq. (3). The trial 
wavefunction is obtained from this by writing 



lO=«MFT>, 

where the Gutzwiller projection operator, 

V G = Y[ n r( 2 - n r), 



(64) 



(65) 



with n r = frafrai enforces the single occupancy 
constraint. Properties of these wavefunctions can be 
calculated numerically using a standard Monte Carlo 
technique. 52 

In Ref. 26, we found that IV'prj) has the lowest energy 
among a class of spin liquid wavefunctions. Furthermore, 
the energy is very low, less than 1% above an extrapola- 
tion of the exact ground state energy to the thermody- 
namic limit; 19 this is achieved without tuning any con- 
tinuous variational parameters to minimize the energy. 
By continuously deforming IV'mft) to include a small 
VBS ordering of the Hastings type and measuring the re- 
sulting variational energy, we argued that ASL is locally 
stable to this type of ordering. Here, we discuss some 
of the properties of this wavefunction. First, we con- 
struct variational excited states to estimate the velocity 
Vf characteristic of the ASL. Next, we proceed to discuss 
the spin-spin and bond-bond correlation functions of the 
wavefunction, and connect them to our understanding of 
the low-energy effective theory. 

We consider finite systems with periodic boundary con- 
ditions for the physical spin operators, so that 



S{r + Lioi) = S(r + L 2 a 2 ) = S{r). 



(66) 



Such a system has 3LiL 2 sites. This still allows the 
consideration of twisted boundary conditions for the 
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FIG. 8: Plot of the mean-field band structure (solid lines) 
and projected band structure (solid rectangles), obtained as 
discussed in the text. The energy is in units of \ J . The band 
structure is plotted along the line from — M2 to M2 in the 
Brillouin zone. The projected band structure can be fit by 
the mean-field band structure, choosing \ = 0-40 ± 0.04. The 
band indices are shown on the right. 

fcrmions - we consider both periodic and antipcriodic 
fcrmion boundary conditions. A technical point that 
plays a role in some of our analysis is that, depending 
on the system size, it is not always possible to construct 
a projected wavefunction that is fully symmetric under 
the kagome space group. The details of this subtlety 
are given in Appendix E. There we show how to con- 
struct symmetric wavefunctions for L\ = L 2 = 4n. The 
technique of Rcf. 52 is easily extended to calculate with 
these wavefunctions, which are superpositions of three 
projected single Slater determinants, and all the correla- 
tion function results shown here are for these fully sym- 
metric wavefunctions. 



A. Fermi velocity 



particle-hole excitation") with respect to the ground 
state as a function of k, q, and the band indices We 
label the bands by indices 1, 2, . . . , 5 as shown in Fig. 8. 
Note that band 5 is doubly degenerate. 

The spectrum of projected particle-hole excitations can 
be translated into a projected band structure. This is 
simplest if there is a single-particle (or single-hole) ex- 
citation that costs zero energy in the mean-field band 
structure. In this case, if the hole in the particle-hole ex- 
citation has zero energy, we can interpret the energy as 
that of the particle, and vice- versa. This situation can be 
arranged by choosing the boundary condition such that 
there are fermions precisely at the nodal points in mo- 
mentum space. Particles and holes created at the nodal 
points (in bands 3 and 4) have zero energy in mean-field 
theory. We study the projected particle-hole excitations 
for L\ = L 2 = 8 with periodic-periodic boundary condi- 
tions, which ensures that we have nodal fcrmions. (This 
boundary condition, on the other hand, has a higher 
ground state energy than that of periodic-antiperiodic 
boundary conditions, in which case there are no nodal 
fermions.) 

To extract the projected band structure, we fix the 
hole at one node in band 3, then scan the momentum 
and band label of the particle in bands 4 and 5 and 
find the energies of this projected state E(\ip(k, q, i,j)}). 
This is essentially a particle excitation and we plot 
E(\ip(k,q,i,j))) - E GS (E GS is the projected ground 
state energy) in Fig. 8 with respect to the fermi en- 
ergy. Similarly we fix the particle at one node in band-4, 
then scan the momentum and band label of the hole in 
bands 1, 2 and 3. This is a hole excitation and we plot 
Egs ~ E(\ip(k,q,i,j))) instead. 



An important parameter in the algebraic spin liquid 
is the "Fermi velocity" vf- Physically, Vf is defined as 
the ratio of energy and momentum for any low energy 
excitation - because the ASL is Lorentz invariant at low 
energies, all excitations propagate with the same veloc- 
ity. In principle, vf can be directly measured via in- 
elastic neutron scattering, by tracking the leading edge 
of the scattering continuum near the T point and the M 
points, where the spin gap is predicted to vanish. Fur- 
thermore, it appears in the coefficients of specific heat 
and magnetic susceptibility, and in various critical scal- 
ing functions that may be measurable. 

To estimate vf, we need to generate low energy exci- 
tations of the ASL. A natural way to do this is to make a 
particle-hole excitation of the mean-field state, then act 
with the projection operator: 

mk,q,i,j))=V G fl t f qj ^ FT ). (67) 

One can calculate the variational energy of the ex- 
cited state \tp(k,q,i,j)) (later referred to as "projected 



It should first be noted that when both k and q lie 
right on the nodal points in bands 3 and 4, the excita- 
tion en- ergy is very small compared to the bandwidth, 
even if the particle and hole are at two different nodes. 
For example when the particle and hole are at the same 
node, we find the excitation energy relative to the ground 
state to be 0.03(1)J; when the particle and hole are at op- 
posite node, we find the excitation energy to be 0.0T(1)J. 
This is consistent with the mean-field band structure, 
where such excitations cost zero energy. We can fit the 
low energy projected particle-hole excitations in bands 3 
and 4 (i.e. the low-energy excitations) with the mean- 
field band structure by tuning the hopping parameter 
X- (Recall that x = t/J is the mean-field hopping in 
units of J.) We find X = 0.40 ± 0.04 (the mean-field 
value Xmean = 0.221 in Hastings' study 13 ). The param- 
eter x determines the fermi velocity vf = y§^7 where 
a is the lattice spacing and J/ks ~ 200K for herber- 
smithite. Therefore, for herbertsmithite, we estimate 
vf ~ 5.0 x 10 3 m/s. 
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FIG. 9: Contour plots of |Soo(<?,d m in)| f° r ( a ) ^min = 0, (b) 
dmin = v3, (c) d m i n = 2, and (d) d m i n = \fl . The scale 
is shown to the right of each plot. The statistical error is 
roughly 1.5 x 1CF 4 , much less than the q — peak height for 
all values of d m m shown. 




FIG. 10: Real space spin correlations measured at a distance 
of half the system size, for L = 4, 8, 12, 16. The straight line 
is a plot of the function 1/L 4 . 

B. Spin-spin correlation function 

We consider the spin-spin correlation function defined 
by 

S ij (R) = ±(S m -S 0j ), (68) 

and its Fourier transform Sij(q) = Yl R e ~ t9R Sij{R)- 
We are primarily interested in the spin correlations at 
long distances; in order to better understand these, it is 
convenient to define the following object: 

S lJ (q,d min )= Yl e-*- R Sij(R). (69) 

\r(R,i)-r(Q,j)\>d m m 

The distance \r(R, i) — r(0,j)\ above is understood to 
mean the shortest distance between the lattice points 
r(R,i) and r(0,j), accounting for periodic boundary 
conditions. Sij(q,d m i a ) is the Fourier transform of the 



spin correlation function, with all pairs of sites with sep- 
aration below d m ; n removed. 

In Fig. 9 we have plotted |5oo(g, d m ; n )| for L\ = L 2 = 
12 at four different values of d m i n . The q = correlations 
dominate as <i m i n is increased; this trend continues for 
larger values of <i m i n . The same is true for 15*10(9, ^min)| 
(not shown). Furthermore, we find that 5*00(0, d m i n ) > 
and 5io(0, (i m i n ) < for all values of d m i n in the L\ = 
L 2 = 12 system (although, for the very largest values of 
rfmin, the result may be dominated by statistical error). 
This suggests that the long-distance q = correlations 
are dominated by the q = pattern of magnetic order 
shown in Fig. 7. 

One expects that the q — spin correlations decay as 
a power law; the exponent of the power law decay can be 
estimated from the behavior of Sij(R) in real space. In 
particular, for a system size L = L\ = L2, we consider 
the behavior of Sij(R) measured at a distance |Jf?| = 
L/2. The 00 component of Sij(R) is plotted for the two 
inequivalent directions R = (L/2)ai, (i/2)(a 2 — a%), and 
the 01 component is plotted for R = (i/2)ai, (L/2)a,2- 
In all four cases, the data are consistent with 

(70) 

This behavior leads to an estimate A m 1=3 2 for the 
monopole scaling dimension (Fig. 10). 

It is perhaps surprising that the M-point magnetic or- 
ders depicted in Fig. 4 do not appear to contribute signif- 
icantly to the spin correlations of the projected wavefunc- 
tion. In particular, if the estimate A m w 2 is accurate, 
we expect the M-point correlations to decay more slowly, 
dominating at long distances. This follows from the ex- 
pectation A at < 2. 28 We wish to make two comments on 
this result. First, we emphasize again that the projected 
wavefunction may not give correct values for critical ex- 
ponents of the ASL. Second, the M-point correlations in 
the wavefunction may indeed decay more slowly than the 
q = correlations, but with a much smaller coefficient. 

C. Bond-bond correlation function 

The bond-bond correlation function has four inequiva- 
lent components, depending on the relative position and 
orientation of the two bonds involved. These are illus- 
trated in Fig. 11, and the corresponding components of 
the correlation function are 

C\R) = ((S Rfi ■ S RA )(S ^o ■ Soa)) - B 2 (71) 

C 2 (R) = ({S RA ■ S R+ai . ) (S*o,o • 5*o,i)) - B 2 (72) 

C 3 (R) = ((5 R ,o-5 R , 2 )(5o,o-5 ,i))-5 2 (73) 

C\R) = ((S R ,2 • 5 R+a2 , )(5o,o • So.i)) - B 2 (74) 

where B = (5o,o ■ 5o,i) is the average dimerization, and 
the subtraction of B 2 ensures that that the correlation 
function decays to zero at infinity. 



FIG. 11: We illustrate the relative bond positions and orien- 
tations that lead to the four inequivalent components of the 
bond-bond correlation function. C'(R) is given by the expec- 
tation value of the product of the dimerization of the black 
bond, with the dimerization of the gray-shaded bond at posi- 
tion i. For example, the gray shaded bond at position 1 gives 
C 1 (ai), and other values of ^(R) are obtained by translating 
the gray bond by a lattice vector. 



FIG. 13: Plot of the coefficients a(L), for L = 8, 12, 16. The 
value of c\ (8) is not shown, as it is negative. Error bars are on 
the order of the symbol size or smaller, except where shown. 
The straight line is a plot of the function e _2 ' 75 /L 2 ' 25 . 
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FIG. 12: Plot of bond-bond correlation function at a distance 
of half the system size. The straight line is a plot of the 
function L~ 2 ' 25 . 



Because computation of bond-bond correlation is ex- 
pensive for larger systems, we only compute some par- 
ticular values at a distance of half of the system size. 
For L\ = L2 = L we consider either R = (L/2)a\, or 
R — (L/2)a2- Table I lists the five sets of data that we 
obtained. We note that, in the second and third rows 
of Table I, the correlation function changes sign (for a 
relatively small value of L) as L increases. This signals 
that these data sets may be further from the scaling limit 
than those that do change sign. Therefore, in Fig. 12 we 
only plot the first, fourth and fifth sets of data, which are 
consistent with power-law decay. We estimate the power 
law to be 



C\R) oc 



1 



R 



2.25±0.05 ' 



(75) 



Surprisingly, this power-law decay does not seem to 
correspond to the Hastings VBS state. To better under- 
stand the long-distance decay of the correlations, we have 



computed C 1 (i?) for a few values of R near (L/2) (01+02) 
- note that, by reflection symmetry, C 1 (-|, -§) = C x (0, -j). 
We assume the long-distance correlations are dominated 
by the three M-points and the T-point, so that, near 
fl=(L/2)(oi + 02), 



c\R)^j2 c ^ L y 



Qi-R 



(76) 



i=0 



where Qo = 0, and Qi for i = 1,2,3 are the three M- 
point wavevectors. The coefficients Cj(L) can then be 



obtained from the four values C 1 ( 



L L\ 
2 ' 2 • 



1. 



C x (f , f + 1) and C^f - 1, f + 1), and are plotted in 
Fig. 13. The coefficient co(L) is significantly larger than 
the M-point coefficients, and behaves consistently with 
the same power-law decay found for C l (R) in Fig. 12. 

The fact that the dominant bond correlations lie near 
q = 0, rather than at the M-points, is puzzling. Note that 
the spin singlet monopoles Wi, even if they are even un- 
der time reversal (t w = 1), carry crystal momenta at the 
M-points and thus do not contribute to the q = correla- 
tions. One possibility is that the dominant long distance 
bond correlations are still not evident for L\ = L2 = 20. 
If the dominant long distance correlations really are at 
q = 0, then cither there is another operator in the ASL 
leading to these correlations that has been missed so far, 
or the projected wavefunction simply does not give a 
good description of the ASL critical behavior. Since, 
to our knowledge, virtually nothing is known analyti- 
cally about the critical properties of Gutzwiller projected 
wavefunctions in two dimensions, it is impossible at the 
moment to decide among these possibilities. 



VII. DISCUSSION 

Here, we discuss our results in the context of experi- 
ments on herbertsmithite, then conclude with a discus- 
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FIG. 14: Herbertsmithite specific heat in zero and 14 Tesla 
magnetic field (courtesy of J. Helton and Y.S. Lee 1 ) together 
with the mean-field theoretical curve. Note that the theoret- 
ical curve is not a fit to the experimental data - there is no 
tunable parameter. We have assumed J/fcs = 200K and used 
the estimate of vf (equivalently, of the hopping parameter \) 
in Sec. VI A. The specific heat is per mole of formula units. 



sion of open theoretical issues. 

The estimate of vp provided by the projected wave- 
function calculations of Sec. VIA allows us to calculate 
the specific heat C(T) and magnetic susceptibility x(^) 
in mean-field theory, and compare to the experimental 
data. The results are 



X(T) 
C(T) 



74.6N A k% 2 
J 2 



(77) 
(78) 



Here, x(T) is the susceptibility per Cu site, and C(T) 
is the specific heat per mole of formula units. The only 
free parameter is the exchange energy J/ks, which we 
take to be 200K. We note that the 1/N corrections to 
the coefficients of C(T) and x{T) have been calculated 
in Ref. 53. Because of the crude nature of our estimate 
for vf, we have not made use of those results in the above 
formulae. C(T) is plotted in Fig. 14 and agrees well with 



the data for T < 30K. 

For the susceptibility, we take the point of view that it 
is dominated by impurities at low temperature. A bet- 
ter estimate of the intrinsic susceptibility of the kagome 
Heisenberg model is given by the NMR Knight shift. In 
both Refs. 4 and 8, a component of NMR spectrum was 
found to have a shift that decreased at low temperature, 
and was argued to correspond more closely to the in- 
trinsic susceptibility. In the case of Ref. 8, saturation is 
observed at the lowest temperatures. In both cases, the 
slope of the decrease is consistent with that of the cal- 
culated within a factor of two. In order to make 
this comparison, one must convert Knight shift to units 
of magnetic susceptibility - this can be done using the 
quoted hyperfine constants in Refs. 4 and 8. The satura- 
tion observed in Ref. 8 may be due to DM interaction or 
coupling of impurities to the bulk kagome layers. 

We note that the physics of the ASL is expected to 
be relevant at temperatures below the spinon bandwidth 
W < J, which we estimate to be W w 0.25J following 
the analysis of Sec. VI A (we define W to be the difference 
in energy from the Dirac points to the top of the lowest 
empty band). Taking J/ks ~ 200K, the temperature 
W/ks ~ 50K sets the correct scale for the onset of the 
downturn in the Knight shift. 

Perhaps the most striking prediction is the presence of 
gapless spin-triplet excitations at the T and M points in 
the Brillouin zone. If sufficiently large single crystals be- 
come available, this could be tested with inelastic neutron 
scattering. Such neutron experiments would also allow 
detection of the magnetic competing orders at both the 
T and M points; a strong test for the presence of the ASL 
would be a verification of the scaling forms for x"(<7> w ) 
given in Sec. IV. The NMR relaxation rate also provides, 
in principle, a probe of the magnetic competing orders of 
the ASL via the power law behavior 1/T\ cx T v , where 
rj = 2 An — 1, assuming An < A m (see Sec. IV). Such 
power-law behavior has been observed 4 ' 8 , with r\ 0.5 
(which corresponds to Ajv ~ 0.75). However, the power 
law seems to depend on magnetic field, 4 suggesting that 
the observed relaxation rate is dominated by magnetic 
impurities. 

As discussed in Sec. IV, the VBS competing order 
may be probed via inelastic X-ray or neutron scattering 
measurements of the lineshape of an appropriate optical 
phonon. Again, this would require single crystal samples 
to be prepared. 

Since the spinons carry entropy, we expect they will 
contribute to the thermal conductivity k. Measurement 
of k will be important in herbertsmithite; because it dis- 
tinguishes between localized and dclocalized gapless exci- 
tations, it should be particularly helpful in understanding 
the observed spin liquid behavior. It will also be impor- 
tant to develop a theoretical understanding of thermal 
conductivity in the ASL. 

All our predictions for herbertsmithite are tempered 
by the twin complications of impurities and DM inter- 
action. To the extent that magnetic impurities play a 
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dominant role, their effect may be reduced by applying a 
large enough magnetic field to polarize them. For exam- 
ple, measurements of specific heat 1 and susceptibility 7 
suggest that magnetic impurities are polarized for a field 
of H = 14 T at temperature T w 2K. Assuming that 
the coupling of magnetic impurities to the bulk kagomc 
spins is characterized by an energy J; mp <C J, then, for 
temperature in the range Ji mp CTCW, the effect of 
impurities on the intrinsic physics of the kagome layers 
should not be severe. However, in this regime one needs 
to focus on probes where impurity and intrinsic contri- 
butions can be separated. 

An understanding of the effects of DM interaction 
would be greatly aided by a better knowledge of its mag- 
nitude. To this end, it will be important to measure the 
anisotropy of x(T) when single crystal samples become 
available. If it is possible to have a reasonable range of 
temperature where D p ,D z <C T <C W, then DM inter- 
action should not be important within this range. For 
ksT < D p ,D z , we argued in Sec. V that DM interac- 
tion will induce spontaneous breaking of time reversal 
symmetry. It would be interesting to look for this at 
low temperature in herbertsmithite. We also note that, 
based on the local structure of the Cu-O-Cu bonds in 
herbertsmithite, we expect D z < D p , which is consistent 
with the results of Refs. 40 and 41. It may be the case 
that D z is small enough that it can be ignored. 

We conclude by mentioning some of the open theoret- 
ical issues relevant to the present study. It is important 
to develop a better understanding of the effects of im- 
purities on the ASL. The physics of single, nonmagnetic 
impurities have been studied in Ref. 54. What is still 
needed is a treatment of magnetic impurities, and an un- 
derstanding beyond the single impurity level. It would 
also be useful to make a systematic study of slave fermion 
mean-field states including DM interaction. 

It would be helpful to understand the critical proper- 
ties of the ASL better for Nf = 4. For example, it might 
be possible to calculate the exponent tjn by numerical 
simulations of the effective field theory. In this paper, 
we have used calculations in the projected wavefunction 
approach to try to estimate critical exponents - a better 
understanding of the criticality in projected wavefunc- 
tions, and its relation to that of the ASL, would also be 
helpful. 
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APPENDIX A: CONTINUUM LIMIT OF THE 
MEAN-FIELD STATE 

In this appendix, we first solve the mean-field Hamil- 
tonian Eq. (3), and discuss its band structure. Then, 
focusing on the low-energy excitations near the Dirac 
nodes, we take the continuum limit and demonstrate 
explicitly the relationship between the continuum and 
lattice spinon fields. The realization of the microscopic 
symmetries (e.g. space group) in the continuum theory 
depends crucially on these results; this is discussed in 
Appendix B. 

We work with the 6-site unit cell as shown in Fig. 2. 
Note that the 6-site unit cell is used in order to ac- 
commodate the background 7r-fiux per hexagon of the 
U(l) gauge field; however, translation symmetry is not 
broken, and the true, or underlying, unit cell has 3 
sites as usual. Unit cells are labeled by Bravais lat- 
tice vectors R — 1n\a\ + 712(12, where n\ and ni are 
integers, and the primitive vectors are 2a\ = 2x and 
a 2 = (l/2)cc + (\/3/2)y. Note that we use the symbol 
R for lattice vectors of the underlying 3-site kagome unit 
cell, and R for the enlarged 6-site unit cell. Within each 
unit cell, sites are labeled as shown in Fig. 2. Dropping 
the spin index, the Hamiltonian may be written 

Wmft = -tJ2 {4 [4i + 4 2 1 + 4i tfjfc ~ 

R 

J R2 l J R+a.2,0 J R-cn+a.2,4 

+ ft [/. +/- 1+ ft [/- + /- 1 
J R3 lJ Ri J R5> J R4 lJ R5 J R+cufi 1 

+ 4,^3-WJ + H.c.} (Al) 

We shall always take t > 0; in fact, this involves no loss 
of generality (see the discussion at the end of this ap- 
pendix) . 

Defining the Fourier transform by 
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where N c is the number of unit cells, we may go to mo- 
mentum space and write the Hamiltonian as Hmft = 

t 'Lkfh H ( k )nfkj> where 
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Here, we have defined K\ = e 2tk ' ai and K 2 = e lk a2 . 

The primitive vectors of the reciprocal lattice are cho- 
sen to be 6i/2 = tt(1,-1/V3) and b 2 = 2tt(0, 2tt/v / 3). 
(Again, the reciprocal lattice primitive vectors for the 
underlying unit cell are bi and b 2 .) The Brillouin zone 
can be chosen as shown in Fig. 1. The two highest en- 
ergy bands of H(k) arc completely flat with energy 2t. 



The lowest empty and highest filled bands meet at Fermi 
points located at ±Q; the Fermi energy is ep = t(\/3— 1). 

We focus on the low-energy excitations near the Dirac 
nodes, and so confine our attention to the two bands that 
touch at ep. We choose the following basis for eigenvec- 
tors of H(±Q) at the Fermi energy, 



(e+) T = ( -e"W24 e «r/8 e ~i*/24 Q £ -3ni/S ^ ( A4 ) 

(e+f = -L(- e ™/ 24 e- 5 "/8 _ e ^/24 ^gii^i/u e -™/8 ) (A5) 

( e -)T = _L^_ e - l7r /24 o e 5x,/8 _g-«r/24 ^-11^/24 e i,/8 ) (A6) 

( e -)T = J_ ( e m/24 ^-13^/24 _ e -«r/8 _ e «r/24 q g-5™/8 ) ) ( A 7) 



so that H(±Q)ef — epef . We want to write an effective 
Hamiltonian for these states, for small deviations of the 
momentum from the nodal points. This can be done 
using first order perturbation theory, which leads us to 
write down the effective Hamiltonian 

[H ± } ij = (ef?D ± (q)ef, (A8) 

Here, D±(q) is given by H(±Q + q) — H(±Q), keeping 
only terms first order in q. The result is 

if ± (g) = i=[ fe r 1 +^r 2 ], (A9) 

which is nothing but the Hamiltonian for massless Dirac 
fermions in two spatial dimensions. Note that the ve- 
locity of the Dirac fermions is isotropic, i.e. it does not 
depend on direction in /c-space. 

We can use these results to define continuum fermion 
fields. Restoring the spin index a, we write 

tfa,±(<f)~ ( { &J f ±Q+q > i '' X } ■ (A10) 
\ ( e 2 )iJ±Q+q,i,a ) 

These fields obey the continuum, second-quantized Dirac 



Hamiltonian 

J^^la^xT 1 +q V T 2 )lp aa . (All) 

We remark that the sign of the spinon hopping t is 
unimportant for all of our results. From the point of 
view of the projected wavefunction, this can be observed 
by noting that t —* —t under a particle-hole transfor- 
mation of the lattice spinons, f ra — ► («cr 2 ) ct( g/^ /3 . Be- 
cause this is an SU(2) gauge transformation, it leaves 
the wavefunction invariant. To understand this from the 
effective field theory point of view, we consider the lattice 
gauge theory Hamiltonian of Eq. (5), and explicitly keep 
track of its dependence on t by writing H e s — H c s(t). 
Noting that both the electric field and vector poten- 
tial are odd under particle-hole transformation, we have 
W e ff(£) — > H e s(—t). Let O be any combination of spin 
operators (it need not be a local operator in space or 
time). O is invariant under SU(2) gauge transformations, 
and, in particular, is invariant under particle-hole trans- 
formation. Therefore the expectation value of O will be 
identical if calculated using either TC c e(t) or H c s(—t), 
and we conclude there is no physical distinction between 
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these two effective Hamiltonians. For simplicity, there- 
fore, we always choose t > as we have done above. 

APPENDIX B: SYMMETRIES 

Here, we outline the procedure for calculating the re- 
alization of microscopic symmetries for the continuum 
Dirac field The basic idea is to diagonalize the mean- 
field Hamiltonian for a finite system, and extract the 
needed information from properties of the nodal wave- 
functions. The necessary manipulations are easily car- 
ried out with standard symbolic or numerical packages 
for linear algebra computations. 

We work with the 6-site unit cell, and consider a sys- 
tem with periodic boundary conditions in the 2a\ and 
a 2 directions, so that R, R + 2L\ai and R + L 2 a 2 are all 
identified. We label sites by r or, equivalently, by the pair 
(R,i). The number of sites is then N s = QN C = 6LiL 2 , 
and the Hamiltonian is a N s x N s matrix defined by 
(H) rr r = =pl for r and r' nearest neighbors, and zero 
otherwise. The negative (positive) sign is taken for the 
thick (thin) bonds in Fig. 2. L\ and L 2 must be chosen 
so that the nodal wavevectors ±Q arc in fact present in 
the Brillouin zone of the finite size system. 

The spin plays no role in these manipulations, so we 
drop the spin index, and we can then think of the con- 
tinuum Dirac field as a four-component object. We write 

*a P (Q = 0) = 5>: p (r)/ P) (Bl) 

r 

where a = +, — is the nodal index, p = 1, 2 is the index 
in the two-component Dirac space. The <& a are the nodal 
wavefunctions, satisfying HQ a = ep^ a , and are given by 

$a P (R,i) = - (B2) 

Now consider a symmetry operation S, with the fol- 
lowing action on the fermion fields: 

S : f r — > Tr/s(r)> (B3) 

where in the present case we can take the gauge trans- 
formation TT r — ±1. This induces the following action on 
the wavefunctions: 

(S$„)(5(r)) =7r r $ a ( r ). (B4) 

This allows us to define the matrix of the symmetry op- 
eration by 

(S)s{r),r = 71V (B5) 

for all r, with all other elements zero. 

Next, we can express the action of the symmetry on 
the nodal wavefunctions as 

S$ a = c a6 $ 6 , (B6) 



fb' 




FIG. 15: 2-fold rotation axes for the symmetries of the 
hexagon. Equivalently, these can be thought of as mirror 
symmetries in the plane. 

where the coefficients c a b can be explicitly computed by 
taking inner products. Translating this into the action 
on the fermion field, we have: 

S : * Q - c* ab * b , (B7) 

which gives us the desired result. 

APPENDIX C: GROUP THEORY OF THE 
KAGOME LATTICE 

1. Outline 

Here, we work out some details of the group theory and 
representation theory of the kagome lattice space group. 
The goal is to understand the action of the space group 
on objects invariant under translations by 2a\ and 2a 2l 
which is true of the competing orders within the ASL, 
and to that end we define a "reduced" space group. We 
use these results to classify all possible site and bond 
ordering patterns invariant under T 2ai and T 2a2 . 

2. Point group 

The point group of the kagome lattice is Dg, which is 
the symmetry group of a regular hexagon. The group Dq 
has 12 elements. We define the elements of D 6 by oper- 
ations on the hexagon shown in Fig. 15. D 6 is generated 
by Cq — R n /3, which is a counterclockwise rotation by 
7r/3 about the z-axis piercing the center of the hexagon, 
and C a = lZ y , which is a tt rotation about the a-axis as 
shown in Fig. 15 (this is equivalent to a mirror symmetry 
in the plane). The group is completely specified by the 
relations C% = C 2 a = (C 6 C a ) 2 = l. 55 

The 12 elements of Dq are 

D 6 = {1, Cq, Cg, Cf , Cg , Cf , C a , C a i , C;,, CV, C c , C c >}. 

(CI) 
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Rep. 


Ci 


r 3 


r 2 


C 6 


Ca 


C a > 


Conj. Class 


Ci 


C t 


r 3 


r 3 




Ai 


1 


1 


1 


1 


1 


1 


Number elts. 


1 


3 


1 


3 


8 


A 2 


1 


1 


1 


1 


-1 


-1 


Representative 


1 


{l|«i} 


{C||0} 


{Cl\ai} 


{C||0} 


B 1 


1 


-1 


1 


-1 


1 


-1 














B 2 


1 


-1 


1 


-1 


-1 


1 


Conj. Class 


c 6 


C a 




Ca' 




Ei 


2 


-2 


-1 


1 








Number elts. 


8 


6 


6 


6 


6 


E 2 


2 


2 


-1 


-1 








Representative 


{C 6 \0} 


{Ca |0} 


{C a \a 2 } 


{CV|0} 


{C a '|ai} 



TABLE II: Character table for A>. 55 The first column labels TABLE III: Conjugacy classes of G s2 , shown with their sizes 
the representation, and the following columns give the values and a representative element, 
of the character on each conjugacy class. 



Here, C a >, C , etc, are rotations by n about the axes 
shown in Fig. 15. Dq has the 6 conjugacy classes 

d = {i}, c 6 = {c 6 ,cf}, ci = {g 6 2 ,g 6 4 }, ci = {ci}, 

C a = {C a ,C b ,C c }, and C a , = {C a ,, C v , C c ,}. The char- 
acter table is given in Table II. 55 The faithful represen- 
tation Ei is obtained by straightforwardly representing 
Z?6 in terms of 2 x 2 rotation and reflection matrices. A 2 
can be obtained from this representation by taking the 
determinant - it distinguishes rotations from mirror sym- 
metries in the plane. Note that Ei = Bi <g> E\ — B 2 ®Ei. 

3. Space group 

The space group of the kagome lattice is generated by 
translations, and the Dq point group operations preserv- 
ing the center of a particular hexagon. A general element 
of the space group can by represented by the Seitz op- 
erator {-R|£}, where R £ 0(2), and t = niai + n 2 a 2 is 
a lattice translation vector. This operator is defined in 
terms of its action on a lattice point r: 

{R\t}r = Rr + t. (C2) 

The entire space group can be generated by {Ge|0} = 
R*/3, {C a \0} = TZ y and {l|oi} = T 0l . 

We now recall some general facts about space groups. 
Denote the space group by G s . The translation group Gt 
is a normal subgroup (i.e. gG t g~ x = G t for all g £ G s ). 
The factor group G p — G s /G t is the point group. 

4. Reduced space group 

Suppose we are interested in the transformation prop- 
erties of some object invariant under translation by a cer- 
tain number of lattice vectors, say by miai and m 2 a 2 . 
So, acting on this object, we have T™ 1 = T™ 2 = 1. Let 
us suppose mi = m 2 = m, and define the group of trans- 
lations leaving our object invariant by 

G tm = j{l|ra(niai + n 2 a 2 )}; ni, n 2 £ z|. (C3) 

Because this set of vectors is invariant under the kagome 
space group, Gt m forms a normal subgroup of G s . (Note 



this would not be the case if mi ^ m 2 .) Therefore we 
can define the reduced space group as the factor group 
G sm — G s /G tm . The reduced space group completely 
describes the action of the space group on our object of 
interest. 

The enhanced fermion bilinears in our kagome alge- 
braic spin liquid are invariant under translation by 2cti 
and 2a 2 . So, to understand their transformations under 
the space group, we consider m = 2 and G s2 . G s2 can be 
interpreted as the kagome space group, acting on objects 
with crystal momentum lying at the V point (q = 0), or 
one of the three M-points. Alternatively, G s2 gives the 
action of the space group on any ordering pattern on the 
kagome lattice with the 12-site 2x2 unit cell (i.e. 2 unit 
cells by 2 unit cells) invariant under T 2ai and T 2a2 . 

We shall work out some properties of G s2 and deter- 
mine its character table. Elements of G s2 can be rep- 
resented again by Seitz operators {i?|t}, where now we 
restrict t = 0,ai,a 2 ,eti + a 2 . The Seitz operators multi- 
ply by the usual rules, except that if one obtains a vector 
t violating the restriction above, we add and subtract 
integer multiples of 2eti and 2a 2 so that the restriction 
is satisfied - this is just the usual way to multiply ele- 
ments of a factor group. G s2 has 48 elements. There 
is a translation subgroup of G s2 , which we denote by 
T 2 = {{l|i},t = 0,eti,a 2 ,ai + a 2 }. T 2 is a normal 
subgroup, and clearly T 2 ~ Z 2 x Z 2 . The factor group 
G s2 /T 2 is just the point group D 6 again. G s2 has 10 con- 
jugacy classes; each class is listed in Table III, together 
with its size and a representative element. 

The character table for G s2 is given in Table C 4. The 
first six representations are obtained from those of Dq, 
exploiting the fact that Dq ~ G s2 /T 2 . To work out 
the properties of the remaining four representations, first 
note that the sum of the squares of their dimensions must 
add up to 36. The only possibilities of dimensions con- 
sistent with this are {3, 3, 3, 3} and {5, 3, 1, 1}. Suppose 
the second possibility occurs, and Fi is a 1-dimcnsional 
representation. Now it must be the case that U Fl (t) = s 
for all translations t £ T 2 , where s = ±1. This is be- 
cause one-dimensional representations must be constant 
on conjugacy classes, and because the translations all sat- 
isfy t 2 = 1. Now we cannot have s = 1, because then we 
would have obtained a distinct irreducible representation 
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1 
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1 


1 
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1 


1 


1 
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1 


1 


-1 


-1 


-1 


-1 


Bi 


1 


1 


-1 


-1 


1 


-1 


1 


1 


-1 


-1 


B 2 


1 


1 


-1 


-1 


1 


-1 


-1 


-1 


1 


1 


Ex 


2 


2 


-2 


-2 


-1 


1 














E 2 


2 


2 


2 


2 


-1 


-1 














Fx 


3 


-1 


3 


-1 








1 


-1 


1 


-1 


F 2 


3 


-1 


3 


-1 








-1 


1 


-1 


1 


F 3 


3 


-1 


-3 


1 








1 


-1 


-1 


1 


F 4 


3 


-1 


-3 


1 








-1 


1 


1 


-1 



TABLE IV: Character table of G s2 . 




FIG. 16: Depiction of the three patterns of site order trans- 
forming according to the Fx representation. Each lies within 
the 12-site unit cell used to construct all possible site and bond 
ordering patterns. Each pattern's crystal momentum lies at 
one of the M-points of the Brillouin zone Mi (see Fig. 1). 
In each case the coefficients c r are zero or ±1. The circles 
represent c r — 1, the squares c r — — 1. 

of D 6 . Therefore s = — 1. However, we have 

-1 = U F HT ai+a2 ) = U Fl (T ai )U Fl (T a2 ) = 1, (C4) 

a contradiction. So it must be the case that all Fi are 
3-dimensional. 

We can construct the representation Fi explicitly - it 
is made up of the 3x3 matrices XJ F > (S) describing the 
action of the space group on N\ and Nfc, which was in- 
troduced in Eq. (30). These matrices can be determined 
using the symmetry transformations of Sec. II B, and the 
character of F\ is then easily obtained. Finally, we obtain 
the remaining representations by taking tensor products 
of F\ with the d = 1 irreducible representations. Specifi- 
cally, F 2 = A 2 <g> F x , F 3 = Si <g> Fx and F 4 = B 2 ® Fx. 

5. Site and bond ordering patterns 

We can use the above group-theoretic results to clas- 
sify all possible ordering patterns with the 12-site unit 
cell shown in Fig. 16. Such patterns are invariant under 




FIG. 17: Depiction of the three F A bond orders. The patterns 
are labeled by their crystal momenta (Fig. 1). The coefficients 
Cb are zero or ±1. Crosshatched bonds represent c& = 1, and 
gray- shaded bonds c b — — 1. 

translations by 2a\ or 2a 2 , and thus transform under the 
reduced space group G s2 . 

We first focus on patterns of order that can be vi- 
sualized in terms of a real field residing on the lattice 
sites. We are primarily interested in collinear spin order- 
ing patterns, and, in this case, the real field should be 
associated with (S*). We define a 12-dimensional vector 
space, where basis vectors are labeled by |r), and r is one 
of the 12 sites in the unit cell. The action of the space 
group on this vector space is given by S\r) — \S(r)}, and 
we have thus constructed a (reducible) representation of 
the space group, which we shall call V s . The matrices for 
this representation can be explicitly constructed, and the 
character is as shown in Table V. The decomposition of 
Vs into irreducible representations is 

V s = Ax®E 2 ®Fx®F 3 ® F A . (C5) 

An ordering pattern of the type we consider here is 
associated with a real linear combination of the |r), 

| site order) = ^c r |r), (C6) 

r 

where the coefficients c r <~ (S*) are the values of the 
ordering field. For each irreducible representation in 
Eq. (C5), if we find a basis of vectors in the form Eq. (C6), 
then we have found the site ordering pattern transform- 
ing in that representation. We are particularly interested 
in the site ordering pattern transforming as the Fx repre- 
sentation, because this transforms identically to some of 
the competing orders in the Hastings ASL (the N\ and 
A^ fermion bilinears). The patterns transforming in this 
representation are depicted in Fig. 16. 

We now carry out the same analysis as above, but for 
ordering patterns that can be visualized in terms of a 
real field residing on nearest-neighbor bonds. There are 
24 bonds in the unit cell - to each of these we associate 
a vector |6), and bond ordering patterns correspond to 
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TABLE V: Characters of V s and 14 representations of G S 2- 




FIG. 18: Depiction of the three bond orders. The patterns 
are labeled by their crystal momenta (Fig. 1). The coefficients 
Cb are zero or ±1. Crosshatched bonds represent Cb = 1, and 
gray-shaded bonds Cb = — 1. 




FIG. 19: Patterns of bond order formed by superposing the 
three F\ patterns (a), and the three F B patterns (b). The 
coefficients Cb are zero or ±1. Crosshatched bonds represent 
Cb = 1, and gray-shaded bonds C(, = — 1. These two patterns 
can in turn be superposed with the uniform state (all Cb = 1), 
to form the Hastings state, with different strengths of the 
three inequivalent bonds in the unit cell. 




real linear combinations of the form 

| bond order) =^c 6 |fe). (C7) 

b 

As before, the action of the space group on this vector 
space is given by S\b) = 15(6)). We refer to the resulting 
24-dimensional representation of G s2 as T4; its character 
is given in Table V. The decomposition into irreducible 
representations is 

V b = A 1 ®B 2 ®E 1 ®E 2 ®F 1 S)F 1 S)F 2 S)F 3 ®F 4 ®F 4 . (C8) 

It should be noted that F\ occurs twice; this is tied to 
the fact that the Hastings state has three inequivalent 
bonds in its unit cell, as discussed below. 

The first F\ irreducible representation of bond orders 
(F-^) consists of the patterns shown in Fig. 17, the sec- 
ond (Ff) of the patterns shown in Fig. 18. The Ff 4 
bond orders can be superposed to form the pattern shown 
in Fig. 19a, and the orders to form that shown in 
Fig. 19b. Formally, this superposition is achieved by 
taking linear combinations in the vector space defined 
by Eq. (C7). The Hastings state can be viewed as a 
superposition of the two patterns in Fig. 19, together 
with the uniform state where all bonds have the same 
amplitude. Depending on the relative weights of these 
three states, one constructs Hastings states with differ- 
ent strengths of the three inequivalent bonds in its unit 
cell (see Sec. Ill B). 




FIG. 20: Spin configuration on up-pointing triangles for two 
inequivalent q — ground states. In both triangles n r = y. 
In triangle (a), the chirality vector c = z, while, in triangle 
(b),c = -z. 



APPENDIX D: ORDER PARAMETER FOR q = 
STATE. 



Here, we construct the order parameter for the q = 
magnetically ordered ground state of the kagome lattice, 
in terms of the complex vector n = n r + im. The q = 
state is a coplanar ordering formed by choosing the vec- 
tor sum of the ordered moments on a single up-pointing 
triangle to be zero, and then translating this triangle to 
fill the rest of the lattice. The q — state is thus com- 
pletely specified by the orientations of spins on a sin- 
gle up-pointing triangle. We choose n r to be equal to 
the ordered moment on the "top" site of the up-pointing 
triangle, as shown in Fig. 20. The remaining moments 
are specified by the chirality vector c, where c ■ n r = 0. 
Moving counterclockwise around the triangle, the spin on 
each site is rotated 120° from the previous one about the 
axis given by c. This is illustrated in Fig. 20. We choose 
the imaginary part n.j = c x n r . 

The action of various symmetry operations on the or- 
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loop-2a 




loop-2b 



FIG. 21: We plot the global loops on a torus. To demonstrate 
the kagome lattice problem in text, we plot loop-1 which is 
along ai direction, loop-2a and loop-2b which are along a 2 
direction. Loop-2a and loop-2b are separated by one unit cell 
spacing. 



der parameter n is now easily worked out. We have 



T ai : 


; n - 


-> n 


(Dl) 


Rtt/3 ' 


: n - 




(D2) 


Tl y : 


; n - 


-> n* 


(D3) 


T: 


: n - 


■+ -n*. 


(D4) 



APPENDIX E: FULLY SYMMETRIC 
PROJECTED WAVEFUNCTIONS 

In this section we discuss how to construct a Gutzwiller 
projected wavefunction which is fully symmetric under 
the space group. First, we show that a simple projected 
wavefunction on a finite system may break lattice sym- 
metry. 

Two projected wavefunctions are identical (up to a 
phase) if one can be transformed into another by a gauge 
transformation. Therefore a projected wavefunction is 
completely determined by the hopping magnitudes on all 
the bonds and the fluxes through all closed loops (assum- 
ing there is no pairing), which are both gauge invariant. 
On a torus the fluxes through all loops are reduced into 
the fluxes inside all plaquettes, and two global loops go- 
ing across the two boundaries (for instance, loop-1 and 
loop-2a in Fig. 21). 

Consider the projected wavefunction for the ASL on a 
kagome lattice torus with L 2 odd; that is, an odd number 
of unit cells in the a 2 direction. In Fig. 21 we show two 
global loops, loop-2a and loop-2b, which are separated 
by one unit cell spacing. That is, they are related by 
the translation T 0l .) Accordingly there are L 2 unit cells 
contained in the cylinder between loop-2a and loop-2b. 
Because there is 7r flux through each unit cell, in total 
we have ttL 2 — n mod 2n flux through the cylinder. As 
a result, the fluxes through loop-2a and loop-2b differ by 
tt. This means the translation symmetry along the a\ 
direction is broken in the projected wavefunction. 

To avoid this problem, we focus on lattices with even 
L\ and L 2 . Repeating the analysis above, we conclude 
that translation symmetry along both a\ and a 2 direc- 



tions is preserved. However, there are still two issues that 
need to be addressed. First, we shall see that, in order to 
construct a symmetric wavefunction, we need to ensure 
that the mean-field Hamiltonian does not have fermion 
states lying precisely at the Dirac nodes. Second, even if 
nodal states are not present, the point group symmetry 
may still be broken. Both these issues are resolved by 
noting that one still has the freedom to choose the fluxes 
through the two holes of the torus. The two choices pre- 
serving time-reversal and translation symmetry are and 
7r. For L\ = L 2 = AN, flux corresponds to periodic 
boundary conditions, and 7r flux to antiperiodic bound- 
ary conditions. On the other hand, for L\ = L 2 = AN+2, 
flux corresponds to antiperiodic boundary conditions in 
the ai-direction and periodic boundary conditions in the 
a 2 -direction. ir flux corresponds to periodic boundary 
conditions in the ai-direction, and antiperiodic bound- 
ary conditions in the (^-direction. 

If there are nodal states in the mean-field Hamilto- 
nian, there is a ground state degeneracy arising from the 
different ways to fill these states, and the resulting wave- 
functions transform nontrivially under microscopic sym- 
metries. Furthermore, we find the ground state energy 
is usually increased in such situations. Therefore, we al- 
ways choose the boundary conditions to avoid the nodal 
fermions. We introduce the notation [P, A], for example, 
to denote the state with periodic boundary condition in 
ai direction and anti-periodic boundary condition in a 2 
direction. We find that for L x = L 2 = AN, [P, A], [A, P] 
and [A, A] lack nodal fermions, and [P, P] has nodal 
fermions. For L x = L 2 = (AN + 2), [P,P], [P, A] and 
[A, P] lack nodal fermions, and [A, A] has nodal fermions. 

We find that in both bases, the three states avoid- 
ing nodal fermions transform into one another other un- 
der point group symmetry, thus forming a three dimen- 
sional representation of the point group. Because the 
Dq point group only has one and two dimensional ir- 
reducible representations (see Table II), this three di- 
mensional representation is reducible. In particular, for 
L\ = L 2 = (AN + 2), the three dimensional repre- 
sentation is A 2 E 2 , and therefore it is impossible to 
use it to construct a symmetric wavefunction. But, for 
L\ = L 2 = AN, we find the representation is AiQ)E 2 , and 
we are able to construct a symmetric wavefunction, be- 
cause A\ is the trivial representation. This symmetrized 
wavefunction, which is a linear superposition of [P, A], 
[A, P] and [A, A] states, is the projected wavefunction 
that we used to study the spin and bond correlations in 
Sec. VI B and Sec. VIC. 

The fact that it is impossible to construct a fully sym- 
metric wavefunction for L\ = L 2 = (AN + 2) from the 
[P, A] , [P, P] , . . . wavefunctions does not mean it is im- 
possible to construct a fully symmetric projected wave- 
function at all. For example, one may be able to use 
states with twisted boundary conditions, and obtain a 
linear combination invariant under both time reversal 
and space group symmetry. This possibility remains 
to be studied. Finally, we remark that the procedure 
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described here cannot construct fully symmetric wave- 
functions on the 36-site cluster that has been extensively 



studied using exact diagonalization. 15 ' 18 20 
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